function  out = int_th( k, r, t, s_, t_, r_, theta )
%INT_TH Summary of this function goes here
%   Detailed explanation goes here
out = sin(theta) .* Gn(k-1, t_, r_, s_) .* g(t - t_, r, r_, theta, s_);

end

